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Abstract 

^> ■ 

^ .We report the final results of a study of the i^(3770) meson using a data sample collected with the KEDR detector at the VEPP-4M 
I 'electron-positron collider. The data analysis takes into account interference between the resonant and nonresonant DD production, 
^-iwhere the latter is related to the nonresonant part of the energy-dependent form factor Fo- The vector dominance approach and 
I "several empirical parameterizations have been tried for the nonresonant F^^{s). 
Our results for the mass and total width of i//{3710) are 
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M = 3779.2 t]j toj loi MeV, 
r - 24 9 ^"-^ ^"-^ MeV 

^ ~ ^^-^ -4.0 -0.6 -0.9 ^' 

where the first, second and third uncertainties are statistical, systematic and model, respectively. For the electron partial width two 
possible solutions have been found: 

+79 +17 +13 , 



(1) r,, = i54!-!-!-eV, 

f 72 +24 +90 
-80 -26 -10 ' 



(2) r - 414 +72 +24 +90 Y 
i ee - tit _80 -26 -10 ^ • 



Our statistics are insufficient to prefer one solution to another The solution (2) mitigates the problem of non-DD decays but is 
disfavored by potential models. 

It is shown that taking into account the resonance-continuum interference in the near-threshold region affects resonance param- 
eters, thus the results presented can not be directly compared with the corresponding PDG values obtained ignoring this effect. 



1. Introduction 

The preceding Letter of this volume is devoted to the mea- 
surement of the il/{2S ) meson parameters in the KEDR experi- 
ment performed during energy scans from 3.67 to 3.92 GeV at 
the VEPP-4M e+e" collider In this Letter we describe the ap- 
plication of the developed tools to the measurement of i^(3770) 
parameters omitting details common for i/'(25') and (/'(3770). 

Since the discovery of the i/'(3770), seven experiments con- 
tributed to the determination of its parameters, nevertheless the 
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situation with the mass, total width and electron partial width is 
still not clear. 

The incomplete compilation of results reported on i/'(3770) 
mass is presented in Table [1] It does not include the results 
of Refs. il2Lll3[l with the analysis of the e+e" — > DD data of 
BES ili] and the e+e" DDy data of Belle 115|] in which 
the ^(3770) electron width has been fixed in the fits causing a 
mass bias. In addition, the bin size in Belle data around i^(3770) 
seems too large for a simple center-of-bin fitting. These works 
encouraged us to employ the vector dominance model in the 
analysis ll III . 

The values presented form three partially overlapping clus- 
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ters. The first one with (M) = 3772.5 + 0.4 Me V comes from the 
analyses in which interference between resonant and nonreso- 
nant DD production has been ignored In ad- 

dition, the analyses assumed the simplest shape of nonresonant 
DD-cross section similar to that for point-like pseudoscalars in 
QED. The statistical uncertainty in this case is small (in iQ] the 
influence of i/'(4040) and higher ^^'s included in the analysis in- 
creases the i/'(3770) mass uncertainty). The second cluster of 
B DDK analyses has <M) = 3775.6 ± 2.3 MeV (the 

result of 101 is not included because of its uncertain status). The 
third, highest mass, cluster is formed by the analyses account- 
ing for interference and gives <M> = 3777.3 + 1 .3 MeV. 

As was mentioned in Section 5.2 of the previous Letter, tak- 
ing into account the resonance-continuum interference is es- 
sential for a determination of the (^(3770) parameters. A close 
DD production threshold significantly increases the importance 
of that. A consideration of the interference effects is one of the 
primary goals of this experiment 

If interference is ignored in a fit of the measured DD or 
multihadron cross section, a bias appears in the growing con- 
tinuum contribution that causes a bias in the resonance ampli- 
tude and a shift of the mass value. The signs of these effects 
depend on the relative position of the interference peak and dip. 
The DD cross section at the threshold is fixed at zero, therefore 
the weights of the more distant data points in a fit are larger 
than those of the less distant ones. Evidence for a dip after the 
DD cross section maximum is visible in all published data with 
large enough statistics (see, for example. Fig. 1 of Ref. |0]), 
therefore, the artificial mass shift should be negative (under- 
counted events move the resonance peak to the left). That is 
exactly what we observe analyzing the published mass results. 

If the result on mass of |j4[] is ignored, the i^(3770) mass 
value obtained in B decays does not contradict neither to 3772.5 
nor 3777.3 MeV. The interference of the resonant and nonres- 
onant DD yields also takes place in this case but the relation 
between them can differ from that in e^e^ collisions, besides, 
the interference effect can be partially compensated by subtrac- 
tion of the combinatorial background. Thus, the intermediate 
mass value does not seem surprising. 

Below we briefly describe the theoretical basis of the analy- 
sis performed, enter some details concerning the analysis proce- 
dure and not covered in the preceding Letter, present the results 
on the i/'(3770) parameters and discuss their systematic uncer- 
tainties and model dependence. 

2. Multihadron cross section in the vicinity of i^(3770) 

A few approaches can be employed to determine the reso- 
nanceparameters using a multihadron cross section data. In the 
Ref. la] the fit of the R ratio was performed, in the Ref. llT6ll 
the efficiency-corrected cross section was analyzed. There are 
many different sources of multihadron events such as the tf/{2S) 
and i/'(3770) production, the light quark production etc., thus 



Table 1: Incomplete compilation of results on 1^(3770) mass. 



Analysis 


M^(3770)[MeV] 


Comments 


MARK-I [1] 


3774.1 +3 


e^e"— > hadrons 


DELCO [2J 


3772.1 +2 


e^e"— > hadrons 


MARK-II [3J 


3766.1 ±2 


e^e"— > hadrons 


BeUe [4J 


3778.4 ± 3.0 ± 1.3 


b^d'Wk^ 


KEDR [5] 


3773.5 ±0.9 ±0.6 


e^e"— > hadrons ''^^ 


BES-11 [6] 


3772.4 ± 0.4 ± 0.3 


e^e^^ hadrons 


BES-II [TJ 


3772.0 ±1.9 


e^e~^ hadrons 


Belle [8J 


3776.0 ±5.0 ±4.0 


B D^^'D^^K^ 


BaBar [9] 


3775.5 ±2.4 ±0.5 


B DDK 


BaBar [IJ 


3778.8 ± 1.9 ±0.9 


e+e-^ DDy^''> 


KEDRyjJ 


3778.0 ± 1.6 ±0.7 


e+e"^ hadrons '''■''^ 



(b) 



- omitted in the latest PDG edition 



(c) 



— the result on S(B D°D°K+) is superseded by d 
prehminary results reported at various conferences 



'■''^ — interference between resonant and nonresonant 
DD production is taken into account 



the variation of the net detection efficiency in the whole ex- 
periment range can exceed 20% {13]. The calculation of the 
net efficiency implies knowledge of the resonance parameters 
and accounting for the interference effects, therefore an itera- 
tive analysis is required. In this work we fit the observed mul- 
tihadron cross section not corrected for the detection efficiency 
which allows iterations to be avoided. 

2.1. Observed cross section and D-meson form factor 

In the energy range from slightly below the il/{2S ) peak to 
slightly above the DDn threshold the variation of the light quark 
contribution to R (Ruds) is small, so that the multihadron cross 
section observed in the experiment can be written as 



obs 

sd*d- 

emp 



"d*d- 



RC , emp , 

cr^^ + cr . + 

lids 

nDD ^iiDD " i/r(3770) ^ 



(1) 



cr 



DDn 



'The result of U was obtained solely to check consistency with the previous 
measurements. 



where cr^'^'s are theoretical cross sections, e's are correspond- 
ing detection efficiencies, and cr'°"''''s are terms treated empir- 
ically as described below. The RC superscript means that the 
cross section has been corrected for initial state radiation (ISR) 
effects, yiDD stands for the direct ^{'illQ) decay to light hadrons, 
the other (super/sub)scripts seem self-explanatory, is a 
branching fraction. All detection efficiencies explicitly entering 
Eq. ([T) can be kept energy independent with sufficient accuracy 
for the event selection criteria employed (see Sec. 13. lb . 

The first four terms have no peculiarities in the whole en- 
ergy range of the experiment, while the last four are responsible 
for the excess of the cross section in the tp{'illO) region. 
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The fourth term of Eq. ([T]i corresponding to the Ught quark 
contribution can be scaled as 1 where a relatively small 
parameter 6 is due to the energy dependence of the detection 
efficiency and radiative corrections. Possible variation of R„ds 
can also contribute to S. This term can be easily removed from 
the consideration in the fit of the cross section provided that the 
6 value is known. The DDn cross section can be treated as a 
small correction. We took it into account using the approxi- 
mately known shape and an additional fit parameter. 

Calculations for o-^f2S) '^t described in the preced- 
ing Letter, a small contribution of the J/t// tail was calculated 
similarly to the il/{2S ) one, for the DD production cross section 
(here and below D stands for or D") one has 



crl£.iW) = J Z^-5{W' VT^) cr^-5{w' VT^) 
X Tix, W'^) G{W, W) dW'dx, 



(2) 



where T^ix, s) is the probability to lose a fraction of s in the 
initial state radiation llSll . G(W, W) describes a distribution of 
the total collision energy, which can be assumed to be Gaussian 
with an energy spread ctw- 

For the charged mode (D^D ) the factor zd*d- describing 
the Coulomb interaction between the mesons produced lfl9ll is 
taken according to Sommerfeld-Sakharov 1E0I.I2I. 221: 



na/Bn* 

"■D^D- = . X 0{W-2mD*). (3) 

For the neutral mode {D^D ") there is no such interaction, thus 



(4) 



the step functions 6{W -2mo) are shown explicitly to simplify 
some expressions below. 

The cross section cr^^ can be expressed via the form factor 
Fu and Z)-meson velocity in the c.m.system Po'- 



(ToDm ^^Pl Worn? , Pd = ^l-4m2/W2. (5) 

To determine the parameters of resonances above the DD 
threshold, their amplitudes should be separated in Fq: 

Fd{W) = Yj + Fd"^"^)^ (6) 



where is the phase of the /-th resonance relative to i^^*. 

For the resonance with the partial widths Tee and F^^ and 
the total width r(W), one has a Breit-Wigner amplitude 



eJiVeela^) (r^5(W)//j^) w 

(the vacuum polarization factor is included in Tee)- 

Considering T{M) as a nominal resonance width and intro- 
ducing the sum of the branching fractions to all non-DD modes 
^nDD' ^'^^ obtains the energy-dependent DD partial width 

r 1 - ^^^^^ '"^^^^ ""^^^^^ ' ^^^^ ' " ^""^^ (8) 

ZiyoooiM)doOD«(M)+ZD*D-(M)dD^D-{M) 



in line with the PDG prescriptions (p. 808 of Ref. 123|]). Here 
do+D- and ^/^jo-qo are the Blatt-Weisskopf damping factors for a 
vector resonance l24ll : 



^DD 



Pd - qoRi), 



(9) 



where Rq represents the meson radius and qo is the cm. mo- 
mentum of the meson qo - /3d W/2. The partial width de- 
pendence according to Eq. (jSj corresponds to the approach of 
Ref. ifioll . Its simplified form was used in the experiments |[ll|2, 
[st]. The approach is somewhat different from that employed in 
Refs. ia, lZD by EES which does not lead to noticeable changes 
of the i/'(3770) parameters. 

The DDn cross section entering ([T]i as a small correction 
can be calculated with sufficient accuracy using 



y6 



F — 



3^2 I^DDn r DDtiI ' 



(10) 



Pod. = V(l +mD)VW2)(l -(mo. -mo)2/W2) . 

The quantity F^^-j^^ is treated as a fit parameter. 

2.2. Nonresonant D-meson form factor 

The nonresonant part of the form factor can be written as 



fT{W) 



1 



|i-no(W) 



■fD(W) 



(11) 



with foiW) - 1 for point-like particles. Here Hq is the vacuum 
polarization operator except the contributions of all resonances 
which are written separately in We remind that the full po- 
larization operator is calculated using the total cross section of 
e^e" — > hadrons that already includes all resonances, there- 
fore use of the full operator IT instead of Ho in the nonresonant 
amplitude leads to double counting of the resonances and thus 
incorrect values of the lep tonic widths (see also the discussion 
in Sec. 5.3 of the preceding Letter). 

There are no precise theoretical predictions for F'^'^iW). 
The model-independent result can be obtained using the expan- 
sions of Re F^*(W) and ImF^*(W) at the point W = M with 
the coefficients free in the fit. Our statistics are not sufficient for 
that, thus we have to rely on some model or use a pure empiri- 
cal approach as in Ref. by BaBar also taking into account 
the resonance-continuum interference. 

The most certain prediction of the form factor can be ob- 
tained with an application of the Vector Dominance Model 
(VDM) to charm production. Standard VDM assumes that the 
inclusive cross section e^e^ — > hadrons at low energy is sat- 
urated by the interfering contributions of the limited number 
of vector mesons. A similar assumption can be accepted for 
the inclusive e^e^ — > cc cross section and its exclusive modes 
such as e^e^ — > DD. The VDM-like analysis of the R ratio 
in the energy range of W - 3.7 H- 5 GeV has been performed 
by BES in Ref. (70, where the light quark contribution R„ds was 
calculated using pQCD. The work cited accounts for tp{377Q), 



i//{4Q4Q), i/'(4160) and i/'(4415) resonances but does not account 
for a possible contribution of i/'(25) decays to DD above the 
threshold. Studies of this contribution Refs. Il2l 11311 include 
a theoretical consideration and some analysis of the DD cross 
section measured by BES as well as by BELLE. In this work 
we employ VDM in a simplified form 



(12) 



where Fq is a real constant representing the contributions of 
the (/'(4040) and higher i/r's. The i/'(25) contribution to the 
DD form factor F'j^^^^ was calculated using Eq. (|7| with the 
D^D^ and D^D partial widths defined similarly to Eq. ([8]) 
with a specific value of the effective radius Rq. The value of 
r^^^\M,,y) = r^l'^*(M,,y) + r^fi*(M,,y) at some reference 
point Mref, as well as the constant Fq, should be obtained from 
the data fit (M„/ = 3778 MeV was used). The partial width 
ratio r^o^J/ro^o' presumably close to that of i/'(3770). 

To evaluate the model dependence of the i/'(3770) parame- 
ters we tried a few nonresonant form factor parameterizations, 
which do not assume vector dominance. The most popular em- 
pirical parameterization is probably exponential: 



where qo is the cm. momentum |25|]. It is well motivated far 
above the threshold but has few parameters to describe the low 
energy region. Instead of it we used 



fo 



8q 



(1 + aqq]j + bqq'l^)" 



(n = 0.5,1). 



(13) 



The minus sign is chosen to match the 4/(28) dominance ex- 
pectations. In the case « = 0.5, bg = 0, the nonresonant cross 
section acquires the Blatt-Weisskopf factor (|9]l with Rq - Uq. 
The case « = 1 corresponds to a more rapid form factor fall. Use 
of two parameters Uq and bq allows us to take into account in the 
limited energy range the increase of the DD cross section de- 
scribed by the G(3900) structure in the Ref. Alternatively, 
the dependence on W-nio 



fo 



8w 



1 -t- awiW-lmo) + bwiW-lmo)^ 
and combined dependences 

8qW 



fo 



(1 + Oqw (W-lmo) + bqw q]))" 



(14) 



(15) 



were considered. 

To check validity of the ip{2S) domination hypothesis in 
Eq. (fT2l) the following parameterizations were used: 



fD 



a„, - W \ a,„ - W 



(« =0,1,3), (16) 



where ySo is the D-meson velocity. They are expansions of the 
Breit-Wigner amplitude with the mass «,„ treated as a free pa- 
rameter, the values of n correspond to different assumptions on 
Tiyf) dependence. In case of \^{2S ) dominance the fitted value 
of fl,„ would be close to M^(2s)- 



3. Data analysis 

3.1. Detection ejficiency determination 

To perform a fit of the observed multihadron cross section 
with Eq. ([1]), it is necessary to know six detection efficiencies 
explicitly entering the equation and the detection efficiency e„rfs 
implicitly contained in the term cr™^ related to the continuum 
light quark production. They were determined from Monte 
Carlo simulation. The efficiency e^,£,2) enters Eq. ([TJ in the 
product with the non-DD branching fraction which is 

rather uncertain. That allows one to assume ~ s^(2S)- 

The event selection criteria, which are different for 2004 
and 2006 scans, and the procedure of the detection efficiency 
determination for the \I/{2S ) decay simulation are described in 
detail in the preceding Letter. The tuned version of the BES 
generator 112611 was employed to the obtained \p{2S ) detection 
efficiency in the vicinity of the peak. The same version of the 
generator with parameters optimal for i/^{2S) simulation was 
used to simulate the i//{2S) and J/tfr tails and the continuum 
uds production. To simulate e^e^ — > DD events, DD pairs 
were first generated with the proper angular distribution. De- 
cays of D mesons were simulated using the routine LU2ENT 
of the JETSET 7.4 package ll27ll . The decay tables of JETSET 
were updated according to those of the PDG review ll23ll . 



Table 2: Detection efficiency for the processes of interest and its 
variation in the experiment energy range AW ^ 200 MeV. 



Process 


£2004 


£2006 


As Is, % 




0.75 ± 0.02 


0.84 ± 0.02 


+ 1.0 ±0.3 




0.74 + 0.02 


0.81 ±0.02 


+ 1.0 ±0.3 


<A(25) 


0.63 ±0.01 


0.72 ± 0.01 


-0.1 ±0.1 




0.50 ± 0.02 


0.60 ± 0.02 


-0.2 ±0.1 


uds 


0.55 ± 0.02 


0.69 ± 0.02 


+2.1 ±0.5 



The detection efficiencies for the processes of interest and 
their energy variations are presented in Tabled The systematic 
uncertainties of the efficiencies s jj^ and Suds were estimated by 
variation of JETSET parameters preserving the mean value of 
the charged multiplicity. The systematic uncertainties on sd+d- 
and ££,0750 were found modifying the decay branching fractions 
of D-mesons within uncertainties quoted in the PDG tables. 

3.2. Fitting of data 

The observed multihadron cross section was fitted as a func- 
tion of W with the expression ([T]i using some assumptions about 
the behaviour of the nonresonant form factor i^^*. The details 
on the likelihood calculation can be found in the preceding Let- 
ter. The following additional constraint was applied 



o(W„/) 



(17) 



D°D 



with the reference mass W„/ : 

served cross section maximum. The value r" - 0.776;^q q25 Il28ll 



: 3773 MeV not far from the ob- 



was used. The world average values were also used for the Jj^ 
mass, total and electronic width. The total width of 1/^(25 ) was 
fixed at the value of 296 + 9 keV obtained in the preceding Let- 
ter. The meson radii of Eq. ^ were fixed at 1 fm and 0.75 fm 
for f(3770) and 1^(25), respectively (Refs. SSIll). Since 
the experimental results on the non-DD fraction of i/'(3770) de- 
cays are controversial and theory expects it to be small, 
we performed the fits with S^^^^ = and 0.16 and assigned 
variation of the parameters to the systematic uncertainties. 
The light quark contribution was parameterized as 



uds 



)Ru 



M- ^ 



<(M^(25)), (18) 



where 6^^^ is a radiative correction of about 0.12, R„ds is a 
light quark contribution to the R ratio averaged over the ex- 
periment energy range and cr*^ in a Born level dimuon cross 
section. The values of S^^^ and Buds are constants correspond- 
ing to W = M^(2S)- The parameter 6 was fixed at 0.187 + 0.046 
with the uncertainties dominated by that of the detection effi- 
ciency variation presented in Table |2l The detailed discussion 
can be found below in Sec. 14.41 

A simultaneous fit of three scans has been performed. Each 
scan has its own free parameters (the energy spread crw and 
Ruds) and has other free parameters common for all three scans. 
Among them are the mass M^(2S), the product of the electron 
width and the branching fraction of its decay to hadrons Y^e x 
Shadr for i/'(25); the mass M, the total width T, the electron 
width Tee and the interference phase for The DDn 

contribution was tuned using the free parameter Fj^-j^^. The non- 
resonant form factor has been controlled by either the free pa- 
rameters r^|'^'(M,.f/) (the 1^(25) partial width above the DD 
threshold) and Fq (constant term of the form factor) or by three 
parameters g, a, b defined in Eqs. ( fT3l l. ( fT4b . ( fTSl ) and ( fT6l l. The 
last but not least free parameter was the interference phase (p. 
The total number of free parameters was either 15 or 16. 

The parameters controlling the nonresonant form factor be- 
haviour have strongly correlated asymmetric statistical errors. 
Instead of them we present below the value of the nonresonant 
DD cross section at the resonance peak cr^^iM) and its error 
obtained in fits with modified sets of free parameters (e.g., the 
{Fq, r^^'^*) pair was replaced with the (Fq, o^^) one). 

The observed multihadron cross section for the scans is pre- 
sented in Fig.[T] The curve represents the vector dominance fit. 
The resulting values of ^{2S ) parameters agree very well with 
those obtained fitting the narrow energy range around tf/(2S) 
(previous Letter). The difference in the mass values is 2 keV, 
the variation of the Tee x Si, product is about 0.3%. As a con- 
sistency check, we estimate R„ds for the three scans. The fitted 
values are 2.33±0.10, 2.25±0.09 and 2.31 ±0.06. The weighted 
average ^„rf, = 2.300+0.046+0.108 OtVA^DoF = 0.49/2) agrees 
well with a similar value 2.262 ± 0.122 published by BES in 
Ref. ll32ll and does not contradict to the result of the BES mea- 
surement Q: R^2.14 + 0.01 + 0.07 atW ^ 3.65 GeV. 

The excess of the multihadron cross section in the (^(3770) 
region is shown in Fig. |2] To calculate the excess, the terms 




3700 



3750 



3800 



3850 

W, MeV 



Figure 1: The observed multihadron cross section as a function 
of the cm. energy for the three scans. The curves are the results 
of the vector dominance fit. The detection efficiencies and the 
energy spreads for the scans differ 

1-4 of Eq. ([T]) obtained by the vector dominance fit were sub- 
tracted from the measured cross section at each point, the resid- 
uals were corrected for the detection efficiency calculated by 
weighting the fit terms 5-8. These terms of the fits are pre- 
sented with the curves. The ignored-interference fit and the fits 
with the anomalous line shapes from Ref. lll6ll are presented for 
comparison. 

3.3. On ambiguity of resonance parameters 

It is known that for two interfering resonances the ambigu- 
ity can appear in the resonance amplitudes and the interference 
phase. A detailed study of that issue can be found in Ref. 0411 . 
In the case of two resonances with constant widths complete 
degeneration occurs: one obtains the identical cross sections 
for two combinations of the amplitudes and phase at the same 
values of the mass and width. 

For the energy-dependent widths there is no complete de- 
generation, however, the likelihood function has local maxima 
on the amplitude-phase plane at slightly different mass and width 
values. A similar situation occurs when a resonance interferes 
with a varying continuum. 

In our case the typical difference in equivalent values of 
the two local minima is very small, -2Aln(L) ^ 0.02, thus a 
certain solution can not be chosen. The variation of mass and 
width for possible solutions is small and neglected below. 
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3740 3760 3780 3800 3820 3840 3860 3880 

W, MeV 



Figure 2: Excess of the multihadron cross section in the i^(3770) 
region. The curves show relevant parts of the fits. The error 
bars correspond to the uncertainty of the measured multihadron 
cross section. All data are corrected for the detection efficiency 
which is different in the three scans. See the detailed explana- 
tion in the text. 
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Scan 2 
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Solution 1 (VDM) 
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Figure 3: Excess of the multihadron cross section in the i^(3770) 
region. Solid and short-dashed curves correspond to two VDM 
solutions. Resonant and non-resonant parts are presented sepa- 
rately. 



4. Results of analysis 

4.1. ^(3770) parameters assuming vector dominance 

In Table |3] we compare the i/'(3770) parameters obtained 
under the assumption of ilj{2S ) dominance in the nonresonant 
form factor for two possible solutions with those extracted from 



the ignored-interference fit and the current world average val- 
ues. The small corrections to residual background given below 
in Table |5] of Sec. I4.3l are not applied to results of the fit. The 
continuum DD cross section cr^— is given without the radiative 

DD ° 

correction factor of about 0.75. The values of the mass and the 
electron width for the ignored-interference fit are in good agree- 
ment with the world average ones, while the value of the total 
width deviates from the average one by 1 .5 standard deviation. 
That is probably due to the statistical fluctuation that occurred 
at the three points of the first scan (see Fig.[2]i. 

Taking into account the resonance-continuum interference 
in DD production improves the chi-square of the fits from 91.1/ 
73 to 74.8/71. The phase of the (/^(3770) ampHtude relative 
to the nonresonant form factor is about 171 and 240 degrees 
for the first and second solution, respectively. The nonresonant 
form factor has a negative real part and a small imaginary one. 
At the 1^(3770) peak i^(25) contributes approximately 70% to 
the total value of the nonresonant form factor If the resonance- 
continuum interference is ignored, the total width is not sub- 
stantially affected, however, the mass shift of about -6.0 MeV 
appears as well as dramatical change of the value and error of 
the electron width. The nonresonant DD cross section in this 
case is underestimated as was discussed in the introduction. 

A large splitting of the T^e values is expected in the near- 
threshold region. Let us illustrate that with an example of the 
area method of the r^.^. determination discussed soon after the 
JItj/ discovery |[35|. The electron width is proportional to the 
area under the resonance curve 



(19) 



(the coefficient k is equal to unity for the energy-independent 
total width), therefore the following expression can be obtained 
in absence of radiative corrections for the case when the con- 
tinuum cross section is small compared to the resonant one: 



a RciM) . 
1 + — , / sm I 

3V B,, ^ 



+ — M cos (h X 



37r 



{W-M)^|Y(M)V{W) Rc(W) dW 



(20) 



(W-M)- + r(W)-/4- 



w 



Here a is the fine structure constant, Rq is the continuum con- 
tribution to R, Bee - the e^e" branching fraction and (p is the 
interference phase. The continuum cross section oc (yJRc )' is 
neglected. 

The left part of (l20b corresponds to the area under the mea- 
sured curve {Tee is obtained ignoring the interference), the right 
part has three terms corresponding to the area under the reso- 
nance curve itself (the true Tee), the curve due to the imaginary 
part of the resonance amplitude (it is also proportional to Y^e) 
and the area of the interference wave due to the real part of the 
amplitude. 

Far enough from the threshold, Rc and T{W) are almost 
constant and the integral is suppressed proportionally to F/M. 
However, for a varying Rq and an asymmetric T(W) near the 
threshold, it grows up to 0.02H-0.15 -\/RciM) depending on the 
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Table 3: 1^(3770) fit results for the vector dominance compared to the ignored-interference case. 



Solution 


M, MeV 


r, MeV 


Tee, eV 


(p, degrees 


r^<5', MeV 

DD 




(r^«,nb 

DD 




1 


3779.3!j;^ 




160!™ 


170.7+ 16.7 




-4 8+^ " 

^■°-3.6 


1.83 + 0.96 


35.7 


2 


3779.3!j;^ 


oc T+4.6 
^-'■■-'-4.0 


420!^^ 


239.6 + 8.6 


11 5+"'-5 


_4 9+3.3 


1.71 + 0.86 


35.7 


i.i. 


3773.3 ±0.5 




249!i 








' -0.07 


7.5 


PDG [23J 


3772.92 + 0.35 


27.3+ 1.0 


265 ± 18 













Table 4: ip{311Q) fits results for alternative assumptions on the nonresonant form factor fo- 



Model 


Mass, total width and Pix^) 


Solution 1 (smaller (p) 


Solution 2 (larger 0) 


Equation 


M, MeV 


r, MeV 




4>, degrees 


Tee , eV 


cr^« , nb 

DD 


(f), degrees 


Tee , eV 


cr^l, nh 

DD 


([T3]in=l 


3779. 1+?" 

— l.O 


94 4+5.0 


32.7 


167.6+16.0 


146!« 


1.82 + 0.76 


243.1 + 9.5 


417!^^ 


1.76 + 0.73 


(US n=0.5 


3779.0![;^ 


25.5!^;° 


33.1 


172.2+ 17.3 


172!^« 


1.59 + 0.86 


241.0+ 15.6 


418!^^ 


1.55 + 0.66 


dH 


3779.0!^;! 


24.4!^-i 


32.7 


167.5 + 21.3 


145 


2.09 + 0.87 


243.1 + 9.5 


418!^^ 


2.02 + 0.86 


([I5]in=l 


3779.0!f^ 


24.4!5-i 


32.7 


167.4 + 20.4 


245+68 

^^-^-49 


2.14 + 0.88 


243.0+ 9.6 


422+^^ 


2.07 + 0.86 


O n=0.5 


3779.0![;^ 


25.2!4;2 


33.1 


172.2 + 21.6 


171+68 
^ ' ^ -65 


1.81+0.88 


241.3+ 11.9 


419!^^ 


1.76 + 0.85 


(UHin^O 


3779.6 + 2.0 


25.3 + 6.6 


31.9 


200.4 + 14.7 


137 + 87 


2.20 + 0.93 


230.3 + 33.0 


461+73 


2.47+ 1.37 


(dSlin^l 


3779.6+ 1.9 


25.3 + 6.3 


31.8 


176.1 + 16.6 


154+ 113 


2.14 + 0.91 


239.4+ 14.7 


433 + 74 


1.96 + 0.96 


([I6]in=3 


3779.1 + 1.7 


25.2 + 4.4 


32.9 


126.0+15.8 


139 + 88 


1.89 + 0.90 


282.0+ 16.9 


501 + 89 


2.54 + 0.91 



assumptions about the energy dependence of T and Rc- The 
closeness to the threshold increases the influence of the inter- 
ference effects by an order of magnitude. The coefficient pre- 
ceding COS0 in Eq.(l20t is about 18 keV in the i^(3770) case, 
the fits give Rc(M) ^ 0.3 with a 40^50% statistical uncertainty. 
Together these circumstances make the area method inappli- 
cable to \l/{2>110). A fit of the cross section is obviously not 
so sensitive to taking interference into account, nevertheless a 
splitting of about 260 eV in Table[3]does not seem surprising. 

The resonant and continuum cross sections for the two VDM 
solutions are presented in Fig. [3] The choice of the true solution 
is essential for determination of the non-DD branching fraction 
of 1^(3770). At the cm. energy of 3773 MeV the resonance 
cross section of 3.8![ 4 nb for the first solution and 9.9![ g nb 
for the second one should be compared with the non-DD cross 
section, which is 1.08 + 0.40 + 0.15 nb according to BBS H 
and -0.01 + 0.08!°^,', according to CLEO |i37tl. The branching 
fraction of about 28% for the first solution seems unreasonable, 
however, that can not be considered as a strong argument in 
favor of the second solution until improvement in the non-DD 
cross section accuracy. 

4.2. Model dependence of results 

To evaluate the model dependence of the if/i'illO) parame- 
ters and to check the vahdity of the vector dominance approach, 



the fits were performed with the alternative assumptions about 
the nonresonant form factor /d(W) described in Sec. 12.21 The 
results of the fits are presented in TableU] A few other assump- 
tions were also tried. 

The amplitude-phase ambiguity was found in all cases con- 
sidered. For each fit we assigned the number 1 to the solution 
with a smaller phase value, while the alternative solution got 
the number 2. The electron width for the first solution was al- 
ways smaller than that of the second one and the values for two 
clusters did not overlap. 

The results obtained assuming dependence of the non- 
resonant form factor as in Eq.(fT3]l almost coincide with those 
for W-nio and mixed dependence in Eq. ( fT4l ) and ( fTSl l because 
of the relatively narrow energy range of the experiment. 

The mass parameter a„, of the parameterizations of Eq. ( fT6b 
n - 0, 1, 3 lies between the (//{IS) mass and the DD threshold 
confirming the i//{2S) dominance. Accepting that the (/'(3770) 
parameters corresponding to the vector dominance model are 
the most reliable, we derive the following estimates for the 
model dependence: 6M - '^^1 MeV, = !°-^ MeV for both so- 
lutions and 6Tee = !^^ (!?°) eV, 5cr^« = !°;4 (+0-8) nb for solu- 
tions 1 (2), respectively. The maximum deviation of parameters 
from the VDM results was taken. The definition of the phase (p 
with Eq.(|6ll allows its model-to-model variation, however, the 
difference with VDM exceeds the statistical uncertainty only in 
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the cases ( fTSI ) « = 0, 3 due to a relatively large imaginary part 
of the nonresonant form factor fitted in these cases. 

We also fitted our data with the anomalous line shapes con- 
sidered in the Ref. lfl6ll by BES where a sum of two noninterfer- 
ing Breit-Wigner cross sections and a sum of two destructively 
interfering amplitudes were referred to as Solution 1 and So- 
lution 2, respectively. The parameters of the amplitudes were 
fixed according to Ref. lfl6ll . the two free parameters were intro- 
duced to correct the general normalization and the shift of the 
energy scale. The i/'(3770) scale correction averaged for two 
shapes is 1 .042 + 0.052 at the energy shift of 0.92 + 0.51 MeV 
which demonstrates rather good consistency of KEDR and BES 
data in general. The chi-square probabilities P{x^) are 25.4 
and 30.3% for the solutions 1 and 2, respectively, compared 
to 35.7% for the vector dominance fit. Both shapes provide a 
better description of the data than the single Breit-Wigner am- 
plitude not interfering with the nonresonant one ("i.i". case in 
Table |3]l due to increase of the resonant yield below 3765 MeV. 
In addition, the destructive interference in Solution 2 reduces 
the resonant yield above 3790 MeV but that does not improve 
significantly the general fit quality because of the growth of 
the peculiarity in the 3765^3780 MeV energy region absent in 
our case. Accounting for the resonance-continuum interference 
with a Breit-Wigner resonance amplitude provides the best fit 
of our data although with our statistics we can not exclude the 
shape anomaly reported reported in Ref lfl6ll . It is worth noting 
that interference of the tpi'illQi) structure with the continuum 
DD amplitude should be considered for any shape assumed. 

4. 3. Correction for residual background 

The residual machine background is about 2% of the ob- 
served uds cross section for the scan of 2006 and five time less 
for the scans of 2004 (Sec. 6.3 of the preceding Letter). The es- 
timated numbers of background events are 445 ± 97 and 24 ± 7, 
respectively, whereas the total number of multihadron events 
selected above the DD threshold is 33678. 

To evaluate the impact of the residual background on the re- 
sulting fit parameters, the background admixture was changed 
in a controllable way. To do so, we prepared a few samples of 
background events passing some loose selection criteria but re- 
jected by the multihadron ones. At each data point ; the number 
of multihadron events A^"''' was replaced with N^''' + f ■ N-^, 
where A^*"*^ is the number of events in the background sample 
The fits with the modified number of events show that the vari- 
ations of all fit parameters are proportional to / in the case 
I/I • A^*"' <& N'"''. Selecting the negative / values at which the 
total number of subtracted events matches the expected back- 
ground admixture and taking into account a small detection ef- 
ficiency change, we obtain the corrections for the fit parameters 
presented in Table |5] The systematic uncertainties quoted in- 
clude those of the background admixture estimate and the vari- 
ation of corrections obtained using different background sam- 
ples. 

4.4. Systematic uncertainties 

The main sources of systematic uncertainty in tfr{3770) pa- 
rameters are listed in Table |6] 



Table 5: Correction to fit results compensating the bias due to the 
background admixture. 



Correction 


Solution 1 


Solution 2 


6M, MeV 


-0.06 ± 0.06 


-0.06 ± 0.06 


6r, MeV 


-0.4 ± 0.3 


-0.4 + 0.3 


% 


-3.9 ±2.9 


-1.5 + 1.1 


DD 


+ 1.5 + 0.5 


+ 1.5 ±0.5 


-2004 „ 


-0.5 + 0.3 


-0.5 ± 0.3 


— 2006 


-2.5 ± 1.0 


-2.5 ± 1.0 



When the resonance-continuum interference is taken into 
account, the multihadron cross section becomes rather sensitive 
to the non-DD fraction of (^(3770) decays. It was varied from 
zero to 0.16 as was mentioned in Section [3T2l The variations 
of the 1^(3770) mass and total width were 0.3 and 0.1 MeV, 
respectively. The shift of the electron width was +8.8% for the 
first solution and -2.3% for the second one. 

The uncertainty on the Rq value used to specify the energy- 
dependent width dill, (ig of about 25% (Refs. L30, JJJ) leads to 
these of 0.3 MeV both in the mass and total width. When the 
interference is ignored, the sensitivity to Rq variations reduces 
by a factor of 3. 

The uncertainties due to that of the branching fraction ratio 
for D"D " and D^D are about 0. 1 MeV for the mass and total 
width. Approximately the same uncertainties are obtained be- 
cause of the D meson masses. The estimates were obtained by 
variation of the values within their errors quoted by PDG. 

To estimate uncertainties due to the inaccuracy of the DDn 
cross section treatment at the edge of the energy range of the 
experiment, we used two methods: shrinking of the fit range 
and assumption of the linear dependence on the D-meson cm. 
velocity instead of the cubical one in Eq. (flOl l. The latter cor- 
responds to variation of the eff'ective interaction radius Rq for 
DDn states from zero to infinity. The variations of the mass, to- 
tal width and electron width do not exceed 0.15 MeV, 0.05 MeV 
and 1%, respectively. 

The systematic uncertainties due to the energy dependence 
of the detection efficiencies shown in Table |2] can be neglected 
in all cases except Suds- The latter together with the energy 
dependence of the radiative correction factor and possible R„ds 
variation determine the power in the expression (fTSl l used to pa- 
rameterize the light quark contribution to the multihadron cross 
section. The radiative correction factor 1 + 6'^^ = 1 .125 ± 0.022 
was calculated according to Ref. 1 1811 using the vacuum po- 
larization data compilation by the CMD-2 group reviewed in 
Ref. 1,38.1 . The error quoted includes the uncertainty of the 
detection efficiency dependence on the mass of the hadronic 
system produced via ISR and that of the vacuum polarization 
data. We explicitly considered the J/tfr tail in the cross sec- 
tion ([TJ, thus the correction factor is 14 + 9% less than that used 
in Ref. ll32ll and its variation in the experiment energy range 
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Table 6: Systematic uncertainties on the i^(3770) mass, total 
width and electron partial width. For the latter the uncertainties 
of two solutions are presented where different. The uncertainty 

on the nonresonant DD cross section is also presented. 



Source 


M[MeV] r[MeV] Y,,[%] o-^§[%] 


Theoretical uncertainties and external data 


precision 






+0.0 
-0.5 


+0.0 
-0.2 


+8.8 /+0 

-0 1 -13 


+0 
-12. 


Rq value in T{W) 


0.3 


0.3 


2. 


1.5 




0.1 


0.1 


0.4 


0.8 


D, D masses 


0.06 


0.04 


0.3 


0.5 


DDn cross section 


0.15 


0.05 


1. 


2. 


Detector and accelerator related uncertainties 


Det. efficiency variation 


0.03 


0.04 


2.4 


5. 


Hadronic event selection 


0.3 


0.3 


3. 


5. 


Residual background 


0.06 


0.3 


2.9 


3. 


Luminosity measurement 


0.1 


0.1 


2. 


2. 


Beam energy 


0.03 








Sum in quadrature 


+0.48 
-0.69 


+0.54 
-0.58 


+ 10.5 /+5.7 
-5.7 /-6.1 


+8. 
-14. 



does not reach 0.1%. The precise R measurements at W = 3.07 
and 3.65 MeV l33l do not indicate essential R„ds variation, thus 
we concluded that the uds efficiency variation dominates in the 
uncertainty of the power \-S. Performing the fits with different 
values of 6 we evaluated the uncertainty of the t//{311Q) param- 
eters as 0.03 MeV, 0.04 MeV and 2.4% for the mass, total width 
and electron width, respectively. Compared to that, the energy 
dependence of e^^j gives only a 0.5% bias of the electron width 
and a few keV shifts of the mass and total width. 

The sensitivity of the mass and width to the criteria of the 
multihadron event selection was checked by changing cuts on 
the energy deposited in the calorimeter and conditions on the 
number of tracks. The results were stable within 0.3 MeV. 
The detection efficiency uncertainty due to inaccuracy of the D- 
meson decay ratios [!230 used for the simulation contributes 2% 
to the electron width uncertainty. The dependence on the choice 
of the selection criteria increases it up to 3%. The sensitivity to 
the event selection criteria is partially due to the influence of the 
residual background. We ignore that and treat the background 
correction as an independent uncertainty source which makes 
the uncertainty estimates more conservative. 

Uncertainties due to the luminosity measurement instabil- 
ity are less than 0. 1 MeV for the mass and width. The accuracy 
of the absolute luminosity measurements discussed in the pre- 
ceding Letter contributes less than 2% to the electron width 
uncertainty. The uncertainty on (/'(3770) mass due to the beam 
energy determination does not exceed 30 keV. 



5. Summary 

The parameters of the )/'(3770) meson have been measured 
using the data collected with the KEDR detector at the VEPP- 
4M e^e^ collider Interference of resonant and nonresonant pro- 
duction essential in the near-threshold region has been taken 
into account. 

Our final results on the mass and width of t//(371Q) are: 

M = 3779.2 MeV, 
p _ 24 9 MeV 

^ ~ ^^-^ -4.0 -0.6 -0.9 

The corrections applied to the fit results are listed in Ta- 
ble|5] The third error arises from the model dependence. It was 
estimated comparing the results obtained under the assump- 
tion of vector dominance in the D-meson form factor (quoted 
values) and under a few alternative assumptions which do not 
imply vector dominance. The quoted model errors do not in- 
clude possible deviations of the resonance shape from the Breit- 
Wigner one with usual assumptions about the total width energy 
dependence, which are predicted, e.g., in the coupled-channel 
model H. 

The result on the i/r(3770) mass agrees with that by BaBar 
also taking into account interference (Ref. 03l) and is signifi- 
cantly higher than all results obtained ignoring this effect. The 
mass values obtained studying B-meson decays by BaBar J^t] 
and Belle |8] are lower but do not contradict to our measure- 
ment. 

We got two possible solutions for the i/'(3770) electron par- 
tial width and the radiatively corrected nonresonant DD cross 
section at the mass of i/'(3770): 



(1) 
(2) 



154 



+79 - 



58 -9 -25 
+72 +24 +90 



eV, 



DD 



1.4 ±0.7 



i ee - tit _8o -26 -10 ^ ' 



+0.1 +0.3 1.^ 
-0.2 -0.2 

+0. 1 +0.6 , 



^§=1.3 ±0.7 !-nb 



The phase shifts of the i/'(3770) ampUtude relative to the neg- 
ative nonresonant amplitude are 171 ± 17 and 240 ± 9 degrees 
for solutions (1) and (2), respectively. 

Most of potential models support the first solution and can 
barely tolerate the second one. The increase of the (^(3770) 
mass according to the BaBar and KEDR measurements implies 
the decrease of the 25 -ID mixing used in potential models to 
rise the electron width value above 100 eV (Refs. 
4211 and the reviews l43l 14411 ). The correct choice of the true 
solution is extremely important for a determination of the non- 
DD fraction of i^(3770) decays. 

Because of the large uncertainty the solution (1) does not 
contradict formally to the previously published results, which 
do not take the interference effect into account, the solution (2) 
is only two standard deviations higher than the current world 
average. However, the qualitative consideration and numerical 
estimates confirm that the impact of the resonance-continuum 
interference on the resulting electron width value is large, there- 
fore the resonance parameters obtained taking into account in- 
terference can not be directly compared with the corresponding 
values obtained ignoring this effect. 
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